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The ambipolar-diffusion theory of star formation predicts the formation of fragments in 
molecular clouds with mass-to-flux ratios greater than that of the parent-cloud envelope. By 
contrast, scenarios of turbulence-induced fragmentation do not yield such a robust prediction. 
Based on this property, Crutcher et al. (2009) proposed an observational test that could po- 
tentially discriminate between fragmentation theories. However, the analysis applied to the 
data severely restricts the discriminative power of the test: the authors conclude that they can 
only constrain what they refer to as the "idealized" ambipolar-diffusion theory that assumes 
initially straight-parallel magnetic field lines in the parent cloud. We present an original, self- 
consistent analysis of the same data taking into account the nonuniformity of the magnetic 
field in the cloud envelopes, which is suggested by the data themselves, and we discuss impor- 
tant geometrical effects that must be accounted for in using this test. We show quantitatively 
that the quality of current data does not allow for a strong conclusion about any fragmentation 
theory. Given the discriminative potential of the test, we urge for more and better-quality data. 
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1 INTRODUCTION 

The ratio of the mass and magnetic flux of interstellar molecular 
clouds has received well-deserved observational attention in re- 
cent years (e.g., Crutcher 1999; Heiles & Crutcher 2005). For a 
cloud as a whole, the mass-to-flux ratio is an important input to the 
ambipolar-diffusion theory of fragmentation (or core formation) 
in molecular clouds (e.g., see Fiedler & Mouschovias 1992, eq. 
[8]; 1993, eq. [lc] and associated discussion). What the ambipolar- 
diffusion theory predicts is the mass-to-flux ratio of fragments (or 
cores) in molecular clouds and how this quantity evolves in time 
from typical densities ~ 10 3 cm -3 to densities ~ 10 14 cm -3 
(Tassis & Mouschovias 2007). Observations have been in excel- 
lent quantitative agreement with the theoretical predictions in that 
the mass-to-flux ratio of cores is found to be supercritical by a fac- 
tor 1 - 4 (Crutcher et al. 1994; Crutcher 1999 and correction by 
Shu et al. 1999, pp. 196 - 198; Ciolek & Basu 2000; Troland & 
Crutcher 2008; Falgarone et al. 2008). By contrast, simulations 
of turbulence-driven fragmentation do not find cores with sys- 
tematically greater mass-to-flux ratios than those of their parent 
clouds (e.g., Lunttila et al. 2008). Therefore, the effort by Crutcher 
et al. (2009) (hereinafter CHT) to measure the variation of the 
mass-to-flux ratio from the envelopes to the cores of four molec- 
ular clouds and thereby constrain cloud-fragmentation theories is a 
much needed observational test. 

The effort by CHT to measure the magnetic field in four cloud 



envelopes yielded mostly nondetections, allowing only the place- 
ment of weak upper limits. Also, the data are suggestive of spatial 
variations of the field in the cloud envelopes. This spatial varia- 
tion must be explicitly treated in the data analysis. Instead, CHT 
performed an analysis based on the overly restrictive (and contra- 
dicted by the data) assumption of uniform magnetic field in the en- 
velope, which minimizes the potentially constraining power of their 
observations. CHT attempt to justify their restrictive assumption 
by claiming that they are testing the "idealized ambipolar-diffusion 
model" that assumes initially straight-parallel field lines in the par- 
ent cloud. Thus, if the data and the data analysis in CHT are taken 
at face value, they at best test an input to a theory, not the predic- 
tion of the theory relating to the variation of the mass-to-flux ratio 
from a core to its envelope, given the field strength and its spatial 
variation in the envelope. As we show below, the geometry of the 
field lines in a parent cloud crucially affects the observed varia- 
tion of the mass-to-flux ratio from a core to the envelope while the 
fundamental prediction of the ambipolar-diffusion theory (that the 
mass-to-flux ratio increases from the envelope to the core) remains 
unchanged. 



In this letter, we present a novel analysis of the OH-Zeeman 
data, applicable also to other sets of data that show intrinsic varia- 
tion of the quantity being measured. 
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Table 1. Magnetic Fields and Errors (in fiG) in Four Cloud Envelopes (data 
from CHT). 



Cloud 




f?2 ± (T2 


B3 ± cr 3 


B4 ± 1T4 


L1448CO 


-9 ± 13 


— 11 ± 6 


-7 ±7 


14 ±8 


B217-2 


-13 ±9 


5 ± 6 


6±8 


9 ±13 


L1544 


-3 ±4 


-1 ±4 


22 ±6 


2 ±10 


Bl 


-16 ±6 


0±7 


-3 ±6 


-10 ±5 



2 DATA ANALYSIS 

The CHT data consist of existing OH Zeeman measurements in 
four molecular cloud cores and of four new measurements in the re- 
gion surrounding each of these cores (in the clouds L1448, B217-2, 
L1544, and Bl). For each observation of an envelope's line-of-sight 
magnetic field Bj , CHT quote an associated Gaussian uncertainty 
a-, . These four values for each cloud envelope are shown in Table 
□ col. 2 - 5 (taken from CHT Figs. 2 - 5). 

2.1 The CHT Analysis 

CHT assign a value to the magnetic field strength in each enve- 
lope, which is obtained from a simultaneous least-squares fit over 
the 8 Stokes V spectra (2 spectral lines at each of 4 positions in 
each envelope). The fit gives a single value of the line-of-sight field 
and a single value of its uncertainty in each envelope. The uncer- 
tainty was calculated under the assumption that there is no intrinsic 
spatial variation of the field strength in each cloud envelope and, 
therefore, any spread in the observed Bj values is attributed to ob- 
servational errors. The CHT values for the envelope fields and their 
uncertainties are shown in Table 2, column 2. 

Using this mean field, CHT calculate what they regard as 
the magnetic flux of the envelope, which, combined with the flux 
in the core, is used to obtain the quantity R defined by R = 
(/corcAKorc/Bcorc) / (hav AV e nv/B env ). The quantity I is the 
peak intensity of the spectral line in degrees K, AV is the FWHM 
in km s _1 , and B is the line-of-sight mean field in /iG. A value 
R = 1 would imply that the mass-to-flux ratio does not vary from 
an envelope to a core in the same cloud, while R 1 would im- 
ply a mass-to-flux ratio greater in the core than in the envelope. 
Since most of the CHT measurements of B in cloud envelopes are 
nondetections, the analysis relies sensitively on the treatment and 
propagation of observational uncertainties to obtain limits on the 
derived quantity R. 

As mentioned above, CHT calculate a mean value of B env and 
an uncertainty on this mean under the explicit assumption that the 
magnetic field in the envelope can be described by a unique B env 
value, which their analysis seeks to constrain. However, the mag- 
netic field in the cloud envelope is not known a priori to have a 
unique uniform value. In fact, the data suggest the opposite (e.g., 
observations 2 and 4 in L1448CO differ by more than 3<r; observa- 
tions 1 and 3 in B217-2 differ by more than 2a; observations 1 and 

3 in L 1 544 differ by more than 4a ; observations 1 and 2 in B 1 differ 
by more than 2cr; see Tabl^D- CHT justify this choice by restrict- 
ing their comparison to what they call the "idealized" ambipolar- 
diffusion theory, assuming that the field lines in the molecular cloud 
envelope are straight and parallel. 

If, as the CHT data suggest, the assumption of zero-spread 
Benv is relaxed, the uncertainties CHT calculate are not the relevant 



ones. A simple example will illustrate the point: Consider a cloud 
envelope in which the magnetic field has a distribution of values 
with mean 10 /iG and spread 5 /iG. An observer makes only two 
measurements of the envelope field, each with uncertainty 0.1 /iG. 
The first measurement gives 10 ±0.1 /iG, and the second measure- 
ment gives 14 ± 0.1 /iG (both very likely). Under the CHT assump- 
tion of zero spread, the mean and associated uncertainty are simply 
the average, B mcan = 12 /iG, and the propagated observational 
error, ff mcan = (53 2 =1 u|) 1 ^ 2 /2 = 0.07 /iG. Clearly, however, 
this S m ean differs from its true value by 2 /iG, not by 0.07 /iG. In 
other words, if there is significant spatial variation of B in a cloud 
envelope, the CHT-kind of analysis grossly underestimates the un- 
certainty on the mean. 



2.2 Straight-Parallel Field Lines in Cloud Envelopes? 

A simple inspection of the CHT raw data, taken at face value, re- 
veals that these four clouds do not have straight-parallel field lines 
in their envelopes. But are such clouds expected on the basis of 
theoretical considerations? Straight-parallel field lines in a parent 
cloud is an idealization in some theoretical calculations that ren- 
ders a mathematically complicated multifluid, nonideal MHD sys- 
tem tractable while capturing all the essential physics of the core 
formation and evolution problem. However, it has never been sug- 
gested that in a real cloud, which is an integral part of a dynamic 
ISM, the envelope field lines will be straight and parallel. Distor- 
tions superimposed on the characteristic hour-glass morphology as- 
sociated with the compression of the field lines during gravitational 
core formation are routinely expected. 

Mouschovias & Morton (1985, Fig. 13) had sketched what 
they regarded as a more realistic field geometry in a molecular 
cloud in which there are several (in that case four) magnetically 
connected fragments. That figure is reproduced here as Figure la. 
This configuration can result from relative motion of the fragments 
(labeled A, B, C, and D) within the cloud, due to the cloud's mean 
gravitational field. The motion of a cloud as a whole relative to 
the intercloud medium will also bend the magnetic field lines in 
an almost U-shape, as shown in Figure lb. One can easily visualize 
lines of sight in Figures la and lb (e.g., the line CC') along which a 
measurement would yield £>i OSiCnv > -Bios, core, although the actual 
field strength in the core is greater than that in the envelope as evi- 
denced by the compressed field lines in the core. By contrast, along 
AA', almost the full strength of the core's magnetic field will be 
measured, but only a fraction of the envelope's field strength will be 
detected. Altogether: (1) An idealization in a theoretical calculation 
should not be mistaken for a prediction. (2) Observations that may 
potentially reveal the geometry of the field lines can and should be 
used as input to build a particular model for the observed cloud (as 
done in the case of Bl by Crutcher et al. 1994, and for LI 544 by 
Ciolek & Basu 2000). (3) The geometry of the field lines cannot 
be ignored in analyzing data from observations that measure only 
one component of the magnetic field (e.g., Zeeman observations) if 
the purpose is to test a theory or discriminate between alternative 
theories. The new analysis of the CHT data in § 2.4 accounts for 
the field geometry suggested by the data themselves. 



2.3 Using Data to Calculate a Cloud's Magnetic Flux 

Unlike the CHT analysis, if the data show field reversals, the posi- 
tive and negative values of the measured Bi OSiC nv must not be alge- 
braically averaged (which is what the CHT assumption of a single 
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Figure 1. Schematic diagrams: (a, top) A deformed flux tube that has frag- 
mented along its length in a molecular cloud (from Mouschovias & Morton 
1985). The deformation can be caused by the relative motion of the frag- 
ments, (b, bottom) Deformation of the field lines threading a cloud caused 
by its motion relative to the surrounding medium. The cloud is shown, for 
simplicity, to contain only one fragment (or core), in the neighborhood of 
which the hourglass shape of the field lines had been established during 
core formation but affected by the cloud's motion. The dashed lines in both 
figures represent different lines of sight, whose significance is explained in 
the text. 



magnetic field value in the envelope imposes on the data) and then 
multiplied by the plane of the sky area of the envelope in order 
to obtain its magnetic flux. If the three (perhaps all four) of the 
observed cloud envelopes exhibited true reversals in the field direc- 
tion (but see § 2.4 below), that would imply a bent magnetic flux 
tube threading each cloud. In such a case, only one algebraic sign of 
the magnetic field (the one corresponding to the greatest absolute 
values) should be considered in estimating the magnetic flux of the 
envelope. Figure 2 and its caption clarify this point. 

2.4 A Self-Consistent Analysis of the CHT Data 

Since both the data and theoretical considerations suggest that 
B env exhibits spatial variations, we reanalyze the CHT data prop- 
erly accounting for this effect and thus generalize the relevance of 
the data to realistic clouds (instead of idealized ones with straight- 
parallel field lines). High-quality data analyzed in this manner can 
potentially discriminate between alternative fragmentation theo- 
ries, instead of just providing geometrical input to theories. 




Figure 2. Schematic diagram of a star (e.g., Sun) that has a dipolar mag- 
netic field and is observed (for simplicity) along lines of sight parallel to 
its equatorial plane. To calculate the magnetic flux threading the star by ob- 
serving its surface field B su[ {, only the mean value of B sur f in either the 
northern or in the southern hemisphere should be used. If both values are 
averaged algebraically (as done by CHT), an erroneous flux value of zero 
will be obtained. 

When intrinsic variation of £?i os in a cloud's envelope exists, 
the spread in the observed values is the convolution of the measure- 
ment error and of the intrinsic spread of 73i os . To account for spatial 
variation of 73i os , a likelihood analysis is needed (see Wall & Jenk- 
ins 2003; Lyons 1992; Lee 2004). We assume that the "true" B los 
follows a Gaussian distribution with mean Bq and intrinsic spread 
(To. This distribution is then "sampled" with N measurements Bj, 
each carrying a (Gaussian) error measurement aj . 

At any specific envelope location, there is a probability 
exp [—{B — Bo) 2 /2oq] /y/inao for the magnetic field to have 
a true value B. If the error of measurement at this same location is 
(Tj, then the probability of observing a value Bj of the field, given 
that its true value is B, is exp [—(B — Bj) 2 /2a 2 ^ /y/2waj. How- 
ever, this is not the only way we could get an observed field value 
Bj, since there are many different true values of the field that might 
yield an observation Bj due to measurement errors. To find the to- 
tal probability for a single observation of Bj, we integrate over all 
possible "true" values of the magnetic field at a single location to 
get (the likelihood for a single observation Bj with observational 
uncertainty aj): 



b = 



dB 



exp [—(B - B J ) 2 /(2a 2 )} exp \—{B - B ) 2 /(2o 2 )} 



2no~n 



The likelihood C for TV observations of 73, 



27TCT0 

(1) 

with individual 



uncertainties Oj to come from an intrinsic probability distribution 
with mean Bq and spread er u is the product of the individual like- 
lihoods, C — IljLi h which, after performing the integration in 
equation Q} and some algebraic manipulations, yields (see Venters 
& Pavlidou 2007) 



C (B , (To) = 




exp 



1 N 



3=1 



(Bj - Bp) 2 
°l + ° 2 



(2) 

Any parameters that are not of direct interest (such as the in- 
trinsic spread cro in this case), can be integrated out of the likeli- 
hood. In this way, we can derive the probability distribution of the 
parameter of interest (Bo) while still allowing for all possible val- 
ues in ao, rather than arbitrarily demanding that cro = (as in the 
CHT analysis). The integrated likelihood is called the marginalized 
likelihood, C m \ this probability distribution can then be used to de- 
rive confidence intervals and upper limits where appropriate. The 
(unnormalized) C m for the four clouds is shown in Figures 3a - 3d. 
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Table 2. Magnetic Fields and Their Uncertainties (in ctG), and Upper Limits on the Envelope Field (in fiG) and Ratio R 



Cloud 


Smean i (Tmcan 


Bmax£ ± C£ 


|Bcnv|(sS 2a) 


|fl|«2<r) 


L1448CO 


0±5 


_ 4 +9 


27 


2.0 


B217-2 


+2 ±4 


+2± 7 7 


22 


2.9 


LI 544 


+2 ±3 




29 


5.0 


Bl 


-8 ±3 


-8±1 


20 


1.1 



Co/. 1: The four observed clouds. Col. 2 & 3: Mean field and its uncertainty as given by CHT and by the likelihood data analysis, respectively. Col. 4 & 5: 
Upper limits on the envelope magnetic field and on the ratio R, from the likelihood analysis. 



C m is derived by numerically integrating equation © over cro for 
different values of Bo, and is shown as a solid curve in the four fig- 
ures; the location of the maximum-likelihood estimate for the mean 
Bo is marked by a heavy vertical line in each figure. 

The maximum-likelihood estimates and associated uncertain- 
ties of Bo for the four CHT clouds are shown in Table[2] The uncer- 
tainties are systematically greater than those quoted by CHT. The 
lcr uncertainties are represented by the widths of the dark shaded 
(solid blue) boxes in Figures 3a - 3d. For comparison we show, as 
cross-hatched boxes, the lcr spreads of the values of Bo that CHT 
quote for the same clouds, based on the same data. 

Upper limits for \Ba\ can also be calculated using C m - The 
2a upper limits for the envelope magnetic field are given in Table 
[2] (The 2cr upper limit is that value of |Bo| f° r which a fractional 
area equivalent to the Gaussian 2a (95.4%) is included under the 
marginalized likelihood curve between — | Bo | and |Bq | . Note that 
this is not the maximum-likelihood value of Bo plus two times the 
error. In addition, C m has much longer tails than a gaussian, and 
hence lcr, 2cr, ... values do not scale linearly.) The values of |Bo| 
and — | Bo | which include between them a fractional area of 2cr of 
the marginalized likelihood for the four clouds are marked with 
heavy down arrows in Figures 3a - 3d. 

To correctly propagate uncertainties onto the derived quantity 
R, we do a full Monte-Carlo calculation to derive the probability 
distribution for the values of R as follows. We repeat the follow- 
ing experiment 10 6 times: we draw /core, Ienv, AVcore, AV e nv and 
B C orc from gaussian distributions with mean and spread equal to 
the measurement and uncertainty quoted in CHT; we draw a mean 
value of B cnv from the marginalized likelihood of the previous sec- 
tion; we combine all the "mock observations" of these numbers to 
produce one value of R. We use the 10 6 values of R produced in 
this way to numerically calculate the probability distribution for 
R. We then calculate the 2cr upper limit on \R\ by requiring that 
the fractional integral of this distribution between — R and +R 
be 95.4%. The 2cr upper limits for |B| are given in Table [2] last 
column. These limits are not very strong: \R\ is constrained to be 
smaller than a few, and for no cloud is the upper limit smaller than 
1 — in sharp contrast to the CHT conclusion, that R is in the range 
0.02 - 0.42 in the four observed clouds. 

In our analysis, we relaxed only one of the CHT assumptions 
(that of lack of spatial variation of B cnv , which is not consistent 
with the data). We have retained the implicit assumption of similar 
orientations of the net B cnv and B corc {vectors), because the data 
do not suggest any particular relative orientation of the two vec- 
tors. A more general analysis that would also relax this assumption 
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Figure 3. Marginalized likelihood (not normalized, solid curve) for the 
four cloud envelopes observed by CHT. In each figure: the location of the 
maximum-likelihood value of Bo is marked with a heavy vertical line; the 
width of the dark shaded (solid blue) box indicates the lcr values of Bo ; the 
width of the cross-hatched box indicates the lcr values for Bo according to 
CHT; heavy down arrows mark the 2cr upper limit on | Bo \ ■ 



would increase the uncertainties on R (although not on B Gnv ) and 
would further part from the CHT conclusions. 



3 SUMMARY AND CONCLUSION 

We have presented a self-consistent analysis of recent OH Zeeman 
observations by Crutcher et al. (2009), and have shown how such 
data can be combined in a statistically robust manner to obtain con- 
straints on the mass-to-flux ratio contrast (R) between molecular 
cloud cores and envelopes. Our analysis extends the constraining 
power of such measurements to realistic clouds, beyond the overly 
restrictive assumption of straight-parallel field lines in cloud en- 
velopes, adopted by CHT. We have shown that the CHT data are 
not of good enough quality to constrain the ratio R and thereby test 
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molecular-cloud fragmentation theories: (i) more integration time 
is needed to reduce measurement errors (now most measurements 
of B env yield only upper limits); and (ii) more measurements (in- 
stead of only four) in each envelope are needed to better constrain 
the intrinsic spatial variation of Benv This kind of observations 
coupled with the method of data analysis we presented in this Let- 
ter has great promise and can lead to significant progress in the 
field of ISM physics in general and in understanding the role of 
magnetic fields in molecular-cloud dynamics in particular. 
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